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SYNOPSIS 
of the 

Ph.D. Dissertation 
on 

NONLINEAR OSCILLATIONS OF 
MODERATELY THICK STRUCTURAL ELEMENTS 

ty 

Pashupati Nath Singh 
Department of Mechanical Engineering 
Indian Institute of Technology, Kanpur (India) 

July 1970 

Large amplitude oscillations of straight heani 
and shallow arch, rectangular and circular plates, a 
douhly curved shallow shell, and a spherical cap have 
been studied in this work. Thin and moderately thick 
structural elements have been analysed. Material of these 
elements have been assumed to be isotropic, homogeneous 
and linearly elastic. The elongations and shears are 
assumed to be small but rotations are moderately large. 
Nonlinear strain- displacement relations based on Karma’s 
assumptions have been used. The effects of transverse 
shear and rotatory inertia have been taken into account 
by the so-called method of internal constraints. Either 
the equations of motion have been derived by the use of 
Hamilton’s principle in the way indicated by Reissner or 



ii 


the equations of motion of the elasticity theory have 
been employed. Fixed inertial reference frame has been 
taken and Lagrangian description of motion has been 
employed , 

By assuming space modes and applying Galerkin 
technique the equations of motion have been reduced to 
ordinary differential equations in -which time is the 
independent variable. The resulting equations have been 
combined into one higher order equation. Fourth and 
higher order derivatives of time have been neglected to 
get one second order equation in time. These have negli- 
gible effect on the fundamental period of vibration, 
which has been obtained in the present work.'^ 

The resulting equation has been integrated 
numerically to obtain the period of flexural vibration. 
Relationship between amplitude and period has been 
studied. Results obtained have been compared, for limi- 
ting cases, with' available results and in general the 
agreement is very satisfactory. Effects like dynamic 
buckling, transition from a plate to a shallow shell 
etc. have been studied. 

While the method is general, the study has been 
limited to structures with hinged ends and only the 
fundamental period of vibration has been obtained. 
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chapter I 


INTRODUCTION 


1.1 GEtTERAL 

Beams, plates, and shells are common structural 
elements and the importance of accurate analysis of their 
dvnamj.c behaviour need not be emphasised. It mam very well 
be said that the most Important structural elements of the 
modem airplane, missile, booster, and re-entry vehicle 
comprise of shells, plate$, and beams. l*Jhen these elements 
are sub''ected to severe operational conditions, the ampli- 
tudes need not be small. Since the linear theories are 
valid onl^r for small deflections, it is necessary to employ 
nonlinear theories which allow for finite deflections. 

In elementary analysis the effects of transverse 
shear and rotatory inertia are neglected. These have 
appreciable effects when the members are not very thin and 
also when there are large deformations. Also, the effects 
become more pronounced in higher modes of vibration. 

The expressions "transverse shear" and "rotatory 
inertia", which occur freauently in the present work are 
more or less self-explanatory. The normal to the middle 
surface of the structural element before deformation may 
not remain normal to the middle surface of the elem.ent 
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after deformation. This is ignored in elementary analysis. 
Apart from this, the elements of a structural member 
perform not only a translatory but also a rotatory motion 
during flexural vibrations. The latter is ignored in the 
classical theories. S. Timoshenko is credited with having 
developed the governing equations which includ.e the effects 
of transverse shear and rotatory inertia in the case of a 
beam, althouigh both the effects were discussed as early as 
1859 by M . Bresse*. There are many ways to include these 
effects in the formulation and one such method, the so- 
called method of "internal constraints" will be employed 
here . 

Approximate solutions of nonlinear equations 
become necessary since most of them cannot be integrated 
explicitly in terms of known functions. To be meaningful 
these solutions must be accurate within certain accepta- 
ble limit and must accomm.odate the characteristic features 
of the physical problem. 

In the problems that have been considered, 
elongations and shears are assumed to be small but rota- 
tions are moderately large. This leads to nonlinear 
strain-displacement relations. On thd other hand the 


* Gours de Hecanique Applique's, Mallet-Bachelier, 
Paris, 1359, p. 126. 
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strains are assumed to maintain linear relationships with 
stresses. Thus the problem becomes geometrically nonlinear 
but remains physically (or materially) linear. The study 
of the effect of material nonlinearity is a separate class 
of pioblems in itself and has not been considered here. 

The structural material has also been assumed to be 
isotropic and homogeneous. It may be pointed out that 
in the present age, there is no dearth of materials with 
the above properties. All the structural elements consi- 
dered have been assumed to be of uniform thickness. 

Free flexural vibrations of thin and moderately 
thick structural elements have been studied. Fundamental 
period ( cori?esponding to the fundamental frequency), which 
is of primary concern to structural engineers, has been 
obtained. Sffects of damping and of coupling of modes 
have been neglected.* The structural elements that have 
been analysed are j straight beam and shallow arch, 
rectangular and circular plates, a doublv-curved shallow 
sbell on rectangular plan, and a spherical cap. These 
are the comnonl^/" used elements in engineering design that 
are subiected to bending. Results obtained have been 
compared, in limiting cases, with known solutions and in 


J.D. Ray and C.W. Bert (Ref. 1) have sho\»jn by 
experimental investigation on a beam that these 
effects are negligible. 
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general the agreement Is very satisfactory. By 
computing the values of fundamental period in a few 
cases it has Been amply demonstrated that By an incorrect 
But frequently errployed application of Galerkin technique, 
the valu.es computed come out to Be erroneous. This fact 
for the case of a sandwich plate has Been pointed out 
By Yi-Yuan Yu and J ai-Lue Lai* . 

1.2 LITERATURE SURVEY 

A survey of literature on free vibration problems 
in beams, plates, and shells reveals that either large 
deflection problems neglecting the effects of trans- 
verse shear and rotatory inertia or small deflection 
problems taking these effects into a,ccount have Been 
studie d . 

The linear analysis of free vj.Brations of 
Beams where transverse shear and rotatory inertia 
effects are considered are dealt with extensively in 
literature (Refs. 2, 3, 4, 5, 6, 7, 8, 9, 10, 11 etc.). 
Large amplitude vibrations of beams have Been discussed 
By a number of authors (Refs. 12, 13, 14, 15, 16, 17 etc.). 




"Application of Galerkin's Method to the 
Dynamic Anal'''’sis of Structures", AIAA Journal, 
April 1967, pp. 792-795. 


o 


12 * 

Sringen derived the governing equations for flexural 
vibrations of thin elastic beams when the displacements 
are finite but the strains are smaJ.l. He assumed the 
motion to be extensional which implied the presence of 
a resultant norma.! force acting on the cross-section 
of the beam and obtained very complicated governing 

1 3 

equations. Wagner foimulated a system of governing 
equations for a similar problem using energy techni- 
ques and disregarding the resultant normal force in the 
cross-sections of the beam. Libai^"^ derived a system 
of governing equations of motion for the above problem 
using a convected coordinate system. He made an addi- 
tional assumption that the component of the instantaneous 
velocity of a point on the middle surface of the beam 
could, under certain circumstances, be neglected as 
compared to the component normal to the surface. By 
using a potential function and assuming the extensional 
strains to be very small compared to rotations, he 

obtained much simpler equations than those given by 

15 

Eiu.ngen. Woodall also solved the problem of finite 
amplitude, free, planar oscillations of a thin elastic 
beam assumj.ng the motion to be inextensional and allowing 


The raised numbers indicate references listed 
at the end of the thesis. 
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the existence of a resultant normal force acting on each 
cross-section of the beam. For a siraplv supported beam 
a finite difference, Galerkin, and (regular) perturba- 
tion solutions have been compared, Ray and hert^ have 
investigated both theoretically and experimentally the 
large amplitude vibrations of a flexible beam with 
pinned ends. They have discussed the extent of vali- 

di.ty of various assumptions that are usually made while 

17 

analysing beam vibration problems. Chu and Herrmann 
and Wu and Vinson^® have investigated the beam problem 
as a special case of the plate problem. 

Large amplitude vibration of rectangular 
pio.tes have been studied by a number of authors 

17 

(Refs, 17, 16, 18 etc.). Chu and Herrmann made 
Ivarman’ s type of assumptions for strains and derived 
the governing equations in terms of components of displa- 
cement for thin plates. They used perturbation technique 

to solve the equations. Wu and Winson extended a part 

ts 

of the work by Nash and Modeer , vho had investigated 
thin plate vibrations in the so-called Eerger’ s field 
(they neglected the second invariant of strain). Using 
Hamilton’ s principle in the way indicated by Reissner, 

Wu and Vinson obtained the equations of motion for 
large amplitude oscillations of transversely isotropic 
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plates. In the solution they have included the effect 
of transverse shear but neglected the effect of rotatory 
inertia. Chauhan and Ashvrell^^ studied the vibration of 
a. thin plate as a special case of a shallow shell. They 
investigated small and large amplitude oscilla.tions of 
a doubly curved shell on square plan. They used 
Rayleigh“Ritz method for the linear problem and attempted 
an approximate analysis for the nonlinear problem by 
mahing energy consideration and finally solving the 
nonlrlnear equation by numerical technique . They ha^’e 
neglected the effects of transverse shear 'and rotatory 
inertia. 

Vibration of shallow arches have been studied 
extensively since they are the simplest structural 
elements about which detailed studies may be feasible 
and which ho-ve many common features with shells. 

Snap-through buckling of curved elements such 
as arches and shallow shells has received considerable 
interest in elastic stability. If the applied loading 
causes the curvature to decrease, the element may 
snap-through to a new equilibrium configuration. The 
new equilibrium of buckled configuration is character- 
ised by having a curvature opposite to the undeformed 
shape , The dynamic analog of this problem has been of 
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considerable interest. Studies in this field include 
loadiA?^ uhj-ch is suddenly applied and maintained 
(Refs. 20, 21), periodic vibrations superimposed on a 
static loading (Ref. 2 3), an.d impulsive loading 
(Refs. 24, 25, 26, 27). In these analyses, the spatial 
form of the displacement function is assumed. Approxi- 
mate solutions are then obtained by using variational 
principles. The essential difference in these studies 
is in the choice of the spatial form. 

Static stability of shallow arches has been,_..^ 

po 

studied in a very complete manner by Fung and Kaplan . 

The corresponding dynamic problem was studied by’H'b'ff and 
qO 

Bruce who have analysed various stability problems by 

using, in most cases, a two mode representation for the 
29 

response. Hsu considered the dynamic stability of 
shallow arches against snap-through when subjected to 
impulsive loads. He took all modes in the response into 
account and investigated the effects of various para- 
meters. Later Hsu^*^ and Hsu et. al.^^ investigated the 

stabjlity of shallow arches against snap-through under 

2l P 6 

timewise step loads. Lock and Humphieys used a 
digital and an analog computer respectively to calcula4fe 
the responses of a shallow arch in order to establish 
stability criteria. Lock used two modes while 

■*' Earlier this work had appeared as a report (Ref. 22). 
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Humplireys used six modes in his analysis. 

Vihrations of plates and cylindrical shells 

have been investigated quite extensively whereas results 

on vibrations of spherical shells are less comprehensive, 

Federhofer^^ and Reissner^^ have obtained approximate 

values of fundamental frequencies for the linear axisy- 

mmetric vibrations of a clamped spherical cap by the use 

34 

of Ritz method. Later Reissner , in dealing with 
fundamental frequency of transverse vibrations of shallow 
spherical shells, justified the omission of longitudinal 
Inertia terms from the governing equations of such 
vibrations. His results for the fundamental frequency 
are in good agreement with those obtained more recently 
by Kalnins^^ from the exact solution of the shallow 
shell equations. The effect of thickness- shear defor- 
mation in the linear problem was taken into account by 
Kalnins^^ for a shallow spherical shell and by 
Deresiewicz and Mindlin^"^ and Deresiewicz'^® for circular 

30 

discs. Also, Reissner solved the problem of axisy- 
mmetric vibrations of circular plates of uniform thickness 
and included the effects of transverse shear and rotatory 

inertia. Hovrever, no numerical results were given. 

♦ 

Koplik and Yu investigated vibrations of homogeneous 
and Sandwich spherical caps with clamped edges and 

The method is based on a generalised Hamilton's 
principle and in particular, the associated variational 
equation of motion for nonlinear case that were used 
for plates and shells by Yu (Refs. 41, 42), 
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included the effects of transverse shear and rotatory 
inertia. They have emphasised the transition from a 
circular disk to a shallow shell - the two problems 
that had been treated so far in a more or less isolated 
manner. Humphreys and Bodner , Budiansky and Roth , 
Archer and Lange"^^, Simitses^, Lock et. al."^^, and 
Huang"^® also investigated dynamic buckling of shallow 

spherical shells under certain kinds of loading. 

47 

More recently Grossman, Koplik, and Yu 
solved the geometrically nonlinear vi.bration problem of 
a shallow spherical cap for various bourdary conditions. 
The formulation is general - the effects of transverse 
shear and rotatory inertia are included. Thereafter 
these effects have been neglected and the equations 

AO 

turn out to be the same as those obtained by Connor . 
These simplified equations have been solved and the 
effects of curvature, amplitude, and edge supports on 
the vibration have been studied. Here also the emphasis 
is on transition from a slightly curved plate to a 
shallow shell. This phenomenon is investigated by 
obtaining the ratio of initial rise to thickness for 
which the type of nonlinearity changes from hardening 
to softening. It may be mentioned that nonlinearity 
in the case of flexural vibration of a circular ring 



(as in the case of a shallow shell) is of softening 

type. This was ohserved, both theoretically and experi- 

49 

mentally for a circular ring by Evensen . Subsequently 
he confirmed experimentally (Ref ,50) that nonlinearity 
in the case of thin cylindrical shell is also of soft- 
ening instead of hardening type. Earlier Yu^^ had 
reduced t’-^e problem of nonlinear vibration of cjrLin- 
drical shells to a nonlinear equation in time, which 
contained both quadratic and cubic nonlinear terms 
(the former being explicitly associated with curvature). 
However, by improper choice of displacement, the effect 
of quadratic nonlinearity vras missed and erroneously 
the nonlinearity v;as found to be of hardening type. A 
better form of displacement has been used in Ref .47 to 
retain the effect of quadratic nonlinearity also. 
However, no difference between positive (outwards) and 
negative (inwards) amplitudes of vibration of the sphe- 
rical cap is made in the above work. Consequently, the 
phenomenon of dynamic buckling is not brought out. 
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1.3 0BJ3CT Ai'JD SCOPE 

In the present work, a method has been sugges- 
ted to obtain the natural frequencies for large amplitude 
oscillations of some moderately thick common structural 
elements like beams, plates and shells. The effects of 
transverse shear and rotatory inertia have been included. 
The material is assumed to be linearly elastic but the 
strain" displacement relations are taken to be nonlinear. 
A fixed inertial reference frame and a Lagrangian 
description of motion have been employed. The problems 
of a rectangular plate, a doubly curved shallow shell 
on rectangular plan, and a spherical cap have been 
an.alysed. The results are specialised for straight 
beam, a shallow arch, and a circular disk. Thickness 
of the elements has been taken to be uniform. In all 
cases except that of a spherical cap, the equations of 
motion have been derived using Hamilton’s principle in 
the manner indicated by Reissner. For a spherical cap, 
the equations of motion of the elasticity theory are 
used. However, the final results obtained for this 
case are compared with the results obtained from equa- 
tions derived by a variational approach. The agreement 
for the limiting case is good. 
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In Chapter II , the problems of the above 
structural elements have been formulated. The equations 
of notion have been obtained in terms of components of 
displacement and those of rotation. Assuming approximate 
mode shapes that satisfy the boundary conditions and 
using ■well-knox'jn Galerkin technique, ordinary nonlinear 
differential equations with time as independent variable 
have been derived. These equations have been combined 
into one higher order equation. Neglecting fourth and 
higher order derivatives, which have negligible effect 
on the fundamental period, one second order nonlinear 
equation in time has been obtained. This has been 
integrated numerically. Although the method is general, 
the computations are restricted to obtaining fundamen- 
tal period of elements with hinged ends. 

Chapter III contains the results computed for 
various thicknesses of elements. The results are 
discussed in detail and are compared in limiting cases 
with available results. In general the agreement is 
very satisfactory. 

Chapter IV lists the conclusions derived from 
the results -computed in this work. 

Recommendations for further work are given 


in Chapter V 
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CHAPTER II 


formulation 


2.1 THE CASS OF A RECTANGULAR PLATE 


In this section, the governing equations are 
derived for a rectangular plate (see Fig.l). These are 
then specialised for the case of a straight beam. 
Reissner’s variational theorem (Ref. 52) extended for the 
dynamic case may be written as : 



r/f F.dV - 
V 


rr (p .u + p .V + p .w)ds 

Si 


- fff ^/2iu^ +v2 + w^)dV I dt = 0 (2.1) 

V 


where t^ to tg is an arbitrary interval of time 
^ ~ ^xx* '^xx ^zz* - zz 


<c 


+ T 

. T 

+ X •'X ; 

xy xy 

xz xz ^yz 

= W 

XX 

^ - 



that part of the surface where stress boundary 
conditions are prescribed 

'^x’Py’Pz ~ Components of surface tractions prescribed 
over 
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Equation (2.1) may be reduced to the form (see, for example 
Ref ,53) : 

r'^2 


J 


^1 - 


r^/ff F.dV - ff (p^u + p^v + P2-w)dSj 
V 


+ Iff (u-f'u + V Sy + w £w)dV i dt =0 

! 

7 

Since the interval t^ to t 2 is arbitrary ; 

^ iff!" F.dV - fj (p u + p V + p w)dS V 
i V 3^ X y 2 J 

+ f [ff (u Su + V. <i V + w. ^w)dV = 0 
V 


( 2 . 2 ) 


(2.3) 


In the course of this work finite deformations have 
been assumed. Elongations and shears are small but rota- 
tions are moderately large. Karman’ s type assumptions 
(see, for example Ref. 54) may therefore be made : 

t = u „ 1/2. (w 

XX ,x 


4yy = Py 


Iz V2.(w )‘ 


(2.4) 


= u + V + w .w 
xy ,y ,x ,x ,y 


"^xz = ^,2 ^,x + 


,2 )y ,y jz 
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To include the effects of transverse shear and rotatory 
inertia i.n the governing equations, the so-called raethod 
of '’’internal constraints” will be used here 5 that is, the 
components of elastic displacement must comply with special 
equations of constraint (Ref. 65). The components of 
displacement are then taken in the form : 

u(x 5 y,z,t) = u^(x 5 y,t) + z. (x,y,t) 

v(x,y,z,t) = VQ(x,y,t) + z.|p(x,y,t) (2.5) 

w(x,y,z,t) = w(x,y,t) 


Therefore, 


^-xx 

^o,x 

2. „ 

+ 

1 / 2 . (w Y 


^o,y 

hy 

4- 

1/2. (w 

5 y 

■^ZZ 

= 0 




■N^ 

■’xy 

+ 

trs 

0 

II 


Z 

( + 

Ixz 

= >>: + w „ 

jX 




tyz 






( 2 . 6 ) 


jJ 


The stresses are taken -in the form (see Fig. 2) : 


Rx 12.z.Mx 


XX h 


h" 


yy 


% 12.Z.M, 


* » 

^zz = 0 


(2.7.1) 

(2.7.2) 


(2.7.3) 


(^22 been taken to be zero for the case of 
free vibration. 
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xy 


T 


•X2 


t 


yz 


Nxy 


* z • ^'Ixy 


h 

+ 

h3 

(2.7.4) 

3Qx 

2'h 


(f)2] 

(2.7.5) 

2h 

fi- 


(2.7.6) 


The analysis is restricted to the case when S^~0. Substi- 
tuting eQS.(2.6) and (2.7) in eq.(2.3), integrating 
through the thickness, and carrying out the variation, the 
following equation is obtained : 

r r ^ WQ r.r ^2 _ Nx ^ 

X 


I f r ^ 1/2. (w )2 - + hi) 

J V I 

•^O O — 


Eh Eh 


2 Ny ■ ^ ^’X ^ 

+ (Vo^y + 1/2. (w^y) - ^ + Eh ^ ^ 

Nyv 

+ (u^ „ + v^ ^ + w .^.w S' N , 


'0)7 o,x Gh 


,x--,y^d -xy 


+ (. 


12 .M 


X 


12. V .My^ 


,x 


Eh^ 


Sh'^ 


)SM 


^x 


12. My 12.iJ.Mx 

+ O ^ ^ ^ 

?7 EhO EhO y 


* ( "',7 + e,x ' 


3^) ? M^y 


Gh'^ 


+ ('^ + w 
+ ( ^ + w 


6.Q. 


ix 


jX ’ SGh ^ ^ ^x 

_ iiii) 

,y 5Gh ^ V 
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" ‘-«7,y - '•’xy,x + \ 


+ (“^XjX “ ^xy,7 '*' ■^•^•■’^o^ 

^ ^~‘V?y ~ ^xy,x ■’■ «h.VQ) S'Vq 

+ (“Nx-^,xx - ^^x,x'^5X - ^V^'^sTY " ^y,y-^,y 


Qt.x “ Q, 


2I'T ,^.W - N„,. „.-W 


-7^,x ■^^''xy‘'",xy “ ^'''xy,x*'"5y 


“ ^Vy,y*^,x .hi.w) S' wj 


‘^x^y 


+ j [(\) + (N^y) S Vq + (Nx‘W^x *3. 

\y‘^,y^ + (M j dy 

^ "b 

3. p 

+ r [ (Ny) .svq + 

3 

"’’ + (M^y)r« J dx = 0 

The forces and cou.ples are derived as ; 


57 ■ ^7 


(2,8) 


TiT *“ 

Eh 

["o,X 


( 1 - v^) 

N = 
7 

Eh 1 

_^0 57 '*’ 

( 1 - V^') 1 

II 

N 

^%,y 

+ V 

o,x 


D ( + 
5 X 

a ) 

^57 


^.VQ^y - 

>''^o,x + 1/2 (w^^) 
+ w ^.w „) 


57" J 

(2.9il) 

r )" 

52 :^ 

(2.9.2) 

(2.9.3) 


(2.9.4) 
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= D O , + V. o/ 


M,„ = 




■xy 12 


?y 


^ -rr 


'?y 


_ _ 5Gh . / _• , 

= “TT” ('=v + w 
“X 6 ?x 

^ _ 5Gh , . ^ ,, 

Qy- — W,y) 


(2.9.5) 

(2.9.6) 

(2.9.7) 

(2.9.8) 


Also, from eq.(2.8) the following equations of motion in 
terms of components of rotation and displacement are 


derived 
a h 


f t j “Do/ - S-CLz,^-1 J _ D ( 1+ ) A + 

J J / ^jXX 2 ^?yy 2 ^,xy 6 


o o 


+ w 


jXX 

^,3 y ■ 


) + 12 J dx dy = 0 

f (• d r .D . aOnia . MiLiU 

! i I 2 P,: 

O O V. 


,x 
a b 


( 2 . 10 . 1 ) 


o o ^■ 
5Gh 


jxy 


+■ 


6 


(f.+ w^y) + dy = 0 


(2.10.2) 


a b ^ _ 

I L ■ 


Sh 


u 


Eh 


V, 


Eh 


o o 


.w „.w 


(1-2.^) 2(1- y) °»xy (l_y2) 


Eh 


w „.w 


Eh 


.u^ 


,x ,xx 2(1- v) ,T 2(1+V ) “"0,77 

Eh_ 


. .w .w ^ + °h u_ 

2(1+ V ) jyy o 




dx dy = 0 


(2.10.3) 



2 % 


a b 



o o 


Sh 


V 


Eh 


11. 


Sh 


.w 


( 1 -V^) 2 ( 1 - >•) o,xy (i„ v 2 ) 


.w 


Eh 


w 


Sh 




Eh 


^yy 2 ( 1 - v) 2 ( 1 + v) ‘ 2 ( 1 + v) 


w 

?y 1 


a b 

f rr 

‘ ^ I 


xjr + f h Jg-v^ j dx dy 


0 


(2.10.4) 


Eh 


(ll„ „.w ,„J 


V Eh 


o 0 


UL (1. ,2) 


. (!!„ ^.W + „.W ) 

Oj^ 53^^ o,y ,xx^ 


— — - (u +v ) w 

d+i.' ) o?y o,x^ ''",xy 


ilh 


(1- V 


■' (u- „„.W y"*" Vo ■yv*'W -tr) (ll ^„.W 

^ o,xx 5 X o,yy ,y 2 ( 1 - y) 


jy 


Eh . /- , TT 'I SGh 

•o,xy-",x^ 2(1+V) ^ o,xx*^5y %,yy-^,x^ “ 6 


+V .W , ) - 


•(< + ^ „) 2i5!k-^ XX x^^ ■*■ 

,x '^jy ,xx ,x 


W .(^f ) ) 

,yy d ^ 


Eh 




2Eh 


."W „.%f ,^. 


2(1„/) 'dxx-' ,y^ dyy^'dx^ ^ (.^_^2) 'dx* ,y 


.w 


5Gh „2 


,xy 


V \r + f h w J S' w ^ dx dy = 0 


(2.10.5) 


Further, the boundary conditions as given by eq.(2.8) are: 
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(i) along X = 0 and x = a 


either 

or 

= 0 

11 

o 

X 

il 

o 

o 

11 

xy 

1 

II 

o 

o 

11 

o 

X 

V =0 

o 

11 

( — 1 

0 

xy 

w = 0 





0 


( 2 . 11 ) 


(ii) Similarly along y = 0 and y = b (2.12) 

Equations (2.10) are in the proper form for the applica- 
tion of Galerkin technique. For a plate hinged along all 
the four edges and vi-hrating in the first mode, the 
following expressions for rotation and displacement 

A 

components can be tahen : 

< = An(t) . Cos- ^ Sin 

A = E,(t) Sin iIJS . Cos 4^ 




Uq = M(t) h 


= h(t) h 


w = P(t) h 


a 

Sin 

Cos 

Sin 


2 rr X 
a 

2Trx 

a 

jr X 


Co; 


2 ny 
b 

2 jTy 


(2.13) 


Sin 

Sxn ■"’'™ " 


These satisfy all the boundary conditions for a hinged 
plate__jnamely ; 
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(i) along X = 0 and x -- a 
M = 0 

X 

P = 0 
^ 0=0 

N = 0 
xy 

■w =0 

(ii) similarly along y = 0 and y = b. 

Then, from eqs. (2.10), the following equations are derived 
after integration* : 


A" 

1 



Ai 

4 

®-2 


4 

clo 

c.) 

P 

= 0 

( 2 . 

. 14 . 1 ) 

Bi,’ 

1 

+ 


Bl 

4 

^2 


4 

bg 

P 

= 0 

(2. 

. 14 . 2 ) 

M" 

+ 


M 

4 

^2 

W 

: 4 


p 2 

= 0 

( 2 . 

. 14 . 3 ) 

H" 

4- 


11 

4 

*^2 

M 

4 

CO 

41 

DO 

= 0 

( 2 , 

• 14 . 4 ) 

pu 

4 


CO 

4 

^^2 

P 

4 

®3 

PM 





4 

®4 

PN 

4 



4 

®6 


= 0 

{ 2 . 

. 14 . 5 ) 


From here onwards B/2(l+i^) is substituted 
for G in this section. 
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v/here „ 


(1- 2(1+>>)B^ (l+V)H^ 


2(1- P ) B 


5 JT 

(l+i> )H 


(1_v2)b 2 2(1+V ) (i+v) h2 


2(1- y)B 

5 jr 

( 1+ V )E 


(1- V 2) (i+y ) b2 

p 

2 JT^ 

(1-V )B 


rr^. H 


4(1- v^) 




(2.15) 


(1- i^^)B2 (l+V ) 


2 rr 

(1- V )B 
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do — 


'1 


-ilii — — (1 + JL.) 
4(1-, 2^ B b3^ 


.TT^h2 


32(1- V^) 
5 

12(1+ y ) 


P C 

B‘^ 3- 


(1 + 


-h 

2(1- V^) 
JT H 

2(1- y^)B 

5 JT 


(1 + 


2 ' 


B' 




i2(i+y )H 


5 rr 


12(1+ V )E.H 


In predominantly flexural vibration, in-plane inertia 
Call be neglected (see Ref. 66) and therefore, M” and N" 
are omitted from eqs. (2.14.3) and (2.14.4). Then all 
the five eqs. (2.14) a,re combined into one higher order 
equation with P as the only dependent variable. Neglecting 
fourth and higher derivatives of P, which have negligible 
effect on the fundamental period (see Refs# 2 and 3 ): 


g^ pn p2 + p» + g^ (pt)2^ P g4 P^ + gg P=0 (2.16) 
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where, 

§2 

§3 

§4 

^5 


(ai + \) (3e^ + 3e333_ + 3e^32) 

(a^b^ - agbg) + (a^+ b^) eg - (a3e5 + b3eg) 


6(a^ + b^) (e^ + egj^ + ©432^ 

(a^b^ - agbg) (e^ + + e^Og) 

~ ^2^2^ ®2 *'®'2^3 ^ ^-3^^ ®5 

+ (agbg - a^bg) e^ 


Cgd3 - G gd^ 
“ °2^2 
Cgdg- c^^dg 

Cidi- cgdg 


(2.17)' 


For small amplitude oscillations, terms containing 
products of P and its derivatives can be neglected 
giving, 

gg P” + gg P = 0 (2.18) 
Corresponding non-dimensional period of flexural vibration 
is obtained as 


T 


L 


2 jr 
■/ §5/ §2 


(2.19) 


This is the period of vibration when effects of transverse 
shear and rotatory inertia are included. 
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An erroneous but frequently employed application 
of Galerkin technique is demonstrated below. It is not 
an uncommon practice to combine the equations derived 
after variation in some convenient foirni and then apply 
Galerkin technique. As for example, let there be n equa- 
tions derived by the use of variational principle ; 


If 


fl (F-l^ 


S F^^ J dx dy = 0 


( 2 . 20 ) 


f ( [*n ’V X *'n] V = 0 

Approximate forms of F_ jF can be assumed and the 

JL 5 il 

integration can be carried out. This is the correct 
Galerkin procedure. But if 


h ( ^1, 


f ( F 
n ^ -^1, 


, F^) = 0 


( 2 . 21 ) 


are combined in such a way that the resulting equations 
are not coeff icj.ent equations of ^ ^ , it is not proper 

to apply Galerkin technique in this form. This can be 
shown by carrying out numerical computations. 

From eq,(2.8), the governing equations can also 
be written in the following form : 
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( 2 . 22 . 1 ) 

( 2 . 22 . 2 ) 

(2.22.3) 

(2.22.4) 

(2.22.5) 

( 2 . 22 . 6 ) 

(2.22.7) 

( 2 . 22 . 8 ) 

(2.22.9) 

( 2 . 22 . 10 ) 

( 2 . 22 . 11 ) 


( 2 , 22 . 12 ) 
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-r-rtr 2N,___«W -o-vj. N-v-^Tr'^-TT _ •5r _ 




!:y*'~,xy ■ ^''X5X*'«5X ‘ "'^xy,x*'",y 


+ ’'xy,y''>3!: + ''V,7-'',y + 5x,x + Qy,y ^ = 0 (2.22.13) 


As done in the previous case, in-plane inertia 
terms from eqs . (2 .22,11) and (2.22.12) can he neglected. 
Selecting a function 0 such that ; 


N 


X 





J 


0 

jXX 


(2.23) 


these two equations are identically satisfied. Then from 
eqs. (2.22.1) through (2.22.3), the compatihilitv equation 
for middle plane strains can he written as (see Ref. 57 ): 

0 + Eh ^ “ 0 (2.24) 

From Eqs. (2.22.4) through (2.22.8) Q^, and 

Q are derived as given hy eqs. (2.9.4) through (2,9,8). 
Substituting these in eqs, (2.22,9) and (2.22.10) 


D. 




+ D(i-^) , + p-a t ^' .) , 

2 2 


• ^,xy 


5Gh 

6 


(< +w - 


12 


= 0 


(2.25) 
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D ft + . B + . 




xy 


SGh ( p+ W ) _ = 


6 ■ '■ J-x- 12 

Ecmation (2.22.13) talces the form : 


0 


(2,26) 


5C-h 


^,xx*^,yy ^jyy-^jxx " ^ ^,xy‘^,x7 ^ e 


(<^2- + ^^y + 7 w) -fh w = 0 


(2.27) 


Eliminating o( and ^ from eqs. (2.25), (2.26), and (2,27): 

r D^(l-V ) D.3 h^(3~V ) ^ ^ <i 5DBh(3- V )t7^ 

I 2 ■ 24 *^t2~ 24(1+ V) 


.. 1 


25E^h^ 


5P Eh'" 


— 20 .w 


5Eh 


<^,xx • ’^,yy + ^,yy*^?^ 
5Eh 


,xy' ,xy i2(l+v) 


~ry w-Phw + 


12(1+ V ) 


f SDEh (1- V ) -,4 _ 25B^h^ 

I 24 (1+ V ) 144(1+ v) 


Sf Eh 


4 ^2- 5 ^ f 


144(1+ V ) 


dt^ j 


I w = 0 


(2.28) 


Approximate forms of w and 0 are selected such 
that eq,(2.24) is identically satisfied. Eauation (2,28) 
is satisfied on the average hy the application of Galerkin 
technique. It may however, be mentioned that ail the five 
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boundary conditions cannot be satisfied since only two 
functions are required to be selected in this reduced form 
of equations. Let the following forms be assumed for w 
and 0 for a plate hinged along the edges ; 


w = f(t).h Sin . Sin -2-2^ 

a a 


(2.29) 


0 = f2(t) 



a^(l- v) 


(x^ + y2) 


+ Cos 


2jt 


X 


3 , 


Cos 



(2.30) 


This part of the calculation is limited to a square plate 
of sides a. Tf-^e boundary conditions satisfied by these 
functions are ; 

(i) along X = 0 and x = a, 


u = 0 

o 

= 0 (2.31) 

= 0 


(ii) similarly along y = 0 and y = a (2.32) 
The boundary condition on (5 along x = 0 and x = a however, 
cannot be satisfied. 
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Substituting these expressions for ¥ and 0 in eq.(2.28) 
and multiplying by Sin Sin ^ , integrating over 

the area of the plate, the following equation is obtained 
after suj.table non-dime nsionali sat ion ; 

a^ f" f2 + ag f" + agCfM^f + f + ag f^ = 0 (2.33) 


whe re , 


TI^( 3- V 

1536(1- V^)(l-V ) 


5 n (3- V ) .H- 
768(1- V^) 


jT^(27-5y 


^2 3456(1- 


5 JT^( 23-11 V ) 25 

3456(1- >‘^)(1+V ) 576(1+ V)^ 


J7®(3-V 5 jr'^O-y )H^ 

+ 

768(1- v^)(l-V ) 384(1- 2) 


(2.34) 


5 JT® 


25 JT^ 


1728(1- V^)(l+y )^ 1728(1- V^)(l+V)^ 


jl^(3-V) 5t7®(3_v)^H4 25 .jt^(3-V) 

^ 2304(1- v^)^ 4608(1--))^)^ 4608(1+ y )^( 1- y ) 

The fourth and higher derivatives of time have been 
neglected in derivi.ng eq.(2.33) as has already been 
explained earlier. 
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2.2 T?3 GASS OF A STRAIGHT BEAK 

From the equations of a rectangular plate, the 
equations for a team (see Fig. 3) can be specialised by 
taking variables connected with y direction and the 
Poisson's ratio V to be zero . 

The stress resultants derived from eqs.(2.9) 


are t 


N 


Eh (u^ ^ + 1/2 (w 


M, 


El. 

■X 12 


^ SGh . , X 


(2.35) 


The variational equations as derived from eqs.(2.l0) 
are : 

r (.Ei£_, +ssh 

6 ?x 12 




^ (<< + w „) + ds = 0 


( 2 . 36 . 1 ) 


a 

^ r" *1 >|t 3i|c 

J (-Sh - Eh w w I dx = 0 ( 2 . 36 , 2 ) 

o ^ ’ 


^ jXK' ^ o 


This is to be done before 


2(1+ V ) 


is substituted for G. 




In-pla.Se inertia term has been neglected. 
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(“Sh u . w - Eh u. 
o ,x ,xx 


o,xx 



5Gh j 

6 -iX' 



w 

,xx 


(w „) 


X 



= 0 


(2,.36.3) 


The boundary conditions along x = 0 and x = a specialised 
from eqs. (2.ll) are ; 
eith er or 

<< = 0 ¥i^ = 0 

^0=0 ° (2.37) 

w = 0 = 0 


For a beam hinged at both the ends and vibrating 
in the fundamental mode, the following forms for<a(, u^, and 
w Can be taken : 

°( = (t) Cos -2^ 

1 a 

U = M (t).h Sin (2.38) 

^ a 

W = P (t) . h Sin^^JL. 

a 

Then from eqs. (2.36) the following eqs. are derived : 

m" + M + bg P^= 0 

pn + P^+ Cg P + CgA + c^PM = 0 


(2.39) 
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whe re . 


2 

jr + 




5 JT 


= 


2 ( 1+ y )H 


= 4 JT 


jr% 


= 


9 jr^f# 


32 


5 jr 


Co = 


2 I2(l+y) 


5 Tf 


12H(1+V) 


C4 = - 


O 

JT^H 


(2.40) 


As in the previous case, in-plane inertia is neglected 
and the three eqs, (2.39) are combined into one equation 
with P as the only dependent variable. Retaining terms 
containing derivative's of P only upto the second order : 

g^ P” P^ + gg pt> + p(pt)2 + p3 + P = 0 


(2.41) 
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whe re , 


§1 


3 4 2 

4 JT E 


2 

JT + 




(1 +V)h 2 l 2 (l+v) 


3 4 ^2 

2 JT H 


(2.42) 


= + 5 7T^ 

^4 4 4(l+iJ) 

4 

^ 5 JT 

^5 l2(l+v ) 

For small amplitude oscillations, terms containing 
products of P and its derivatives can be neglected giving, 
gg P" + gg P = 0 (2.43) 

Corresponding non-dimensional period of flexural 
vibration, 


T = 


2 JT 


%/g2 


( 2 ♦44) 
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2.3 THE CASE OF A DOUBLY CURYSD SfiELL ON BSCTAI'IGULAR PLAN 

A shallow shell (see Fig. 4), whose middle surface 

“yc y t * "V" 

is given by WCxjy) = c Sin . Sin is considered 

here, c is the maximum initial rise in unstressed condition. 

For a shell to be shallow c/b should not be greater than 

■i. For such a shell, the first fundamental form of the 
5 

surface can be taken as ds^ = dx^ + dy^ and any integra- 
tion required to be carried over the surface of the shell 
can be carried over the projected plan area without much 
error (Refs. 58, 59), This means that projected forces 
and moments can be taken as the actual forces and moments. 
These are the usual simplifications done in the case of 
shallow shells. 

The components of displacement and distribution 

of stresses through the thickness are assumed as those 

given by eqs. (2.5) and (2.7). 

The components of nonlinear strains are taken on 

SO 

the basis of the derivation given by Marguorre for 

shallow shells, namely, 

= U _ + 1/2 (w ^ TT* -JT 

^xx ,x ' ^ ,x ,x ,x 

^77 = ^,y + 1/2 (w^y)^ '^y* ■^jy 

^ZZ ~ 0 


(2.45) 
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xy 




V 


X 




w 


4- w 


X 






+ vJ 


jy* 


w 


= <=< + W 


X 


Yy2 = ^ W 


?y 


Substituting for u, v, and w : 


o c XT 

= '^n.v + 2 • + 1/8 (W.^) + 


'^0,X 


,x ''5X a 


.w 


X 


Cos -Xj . sin 

a b 


^JY ^o,y 


+ z . /3 + 1/2 (w^y)^ + w 


jy 


Cos ^ . Sin - 21 x 
b a 


^ zz 

- 0 





^xy 

= ^o,y-^ 

’'OjX + 2 • 

( o< + 

jy 

^,x) 

+ W TT* 

jX 


+ c jy 

a 

\y 4 

-S . Sin 

7 T y 
b 

+ w 

b- 


. Sin - 

^ ^ . Cos - 

a 

Jrj 

b 




Yyz = ^ + 


By the variational approach, as in the case of the 
rectangular plate, the force-displacement relations are 
obtained as : 
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N = 


(l-vSy 

> 

+ 

M 

o 

^o,y 

+ cJT . 

■w „ . Cos 

,X 

iLx 

a 

a 

Sin £lx 

. Cos ^ 

1 

a 

b 

i 


v2 . V 


2 


. Sin ^ w ^ . 

b b iY 


ii;n_ r .2 v p 

(1-^^) ^ ^0,x 1/2 (V7^y) + — 


+ . Sin Ee 

0 a 


'TT y + >^0 JT 


a ’■ 


Cos ES. , Sin 7 

3. -j 


^xy= Gh 


[^^0,7 ’'^o,x ■*■ ^,x • ^,y ■*■ • 

Cos ^2 Lx ^ gin . w . Sin iL2^ . Cos iI2l 

a b b jx a b I 


= D ( 


(2.47) 


My = D ( a,y-^ V 


^xy 12 


^ 'Z jY 


_ _ 5Gh / / \ 

^x 6 ^ ^,x^ 


_ — 3Gh f Q , - A 

Qy “ 6 ^ ^ ^,y^ 
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The hounclarv conditions as derived are ; 


(i) along X 

= 0 and X = a 


either 

Qx 


o( = 0 

= 0 


^ = 0 

(1 

o 


O 

li 

O 

- 0 

(2.48) 

o 

1! 

V = 0 


u = 0 

Qx -W* ^%y 

+ U . Cos ^ Sin 

ax a b 

+ ^ N Sin lU. , Cos 41 

b a b 

c: 0 


(ii) similarly along y = 0 and y = b (2.49) 

The components of rotation and displacement, for 
a shell hinged along the four edges, are taken to be the 
same as in eqs. (2.13) for the fundamental mod.e of 
vibration. 

Ultimately the governing eouations in non- 
dimensional form (corresponding to eqs. (2.14) for the 
rectangular plate) reduce to ; 
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A” + ®-2®l 


= 0 


B!^ + + b2% + ^3? = 0 

M’‘ + c^M + C2l'‘T + ^3?^ + C4P = 0 


B'* + d-,B + dgM + doP^ + d.P 


0 


P” + e^P^ + 63 ?^ ® 3 ^ ■*■ ■*■ 


+ e^K + ©gA^ + ®9*®1 “ ^ 


where 


? 



^ (1-V ) 2(1+ V)B^ (!+'>' )H' 


JT' 


2(l-y)B 


a 


5 yr 


■3 (1+v) H 


b. = 


: + 


JT 


t 


^ (l-V^)B^ 2(1 +v) (i+v)h2 


jr 


^2 2(l-v) B 


b^ = 


5 yy~ 


2 (i+y)B.H 


(2.50) 
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4 Jt 2 

c = ^ — -” 5 “ + 7 




2jr 
(1-V )B 


JT^H 


'3 


( 1 + ---) 


4(1- 2 


0. = - Judj. (1 . 1 ) 


2 ( 1 - v^) B" 


d, = 


4 jr" 


2 ^ 


^ (i-v^)b 2 ( 1 +V) 


^ “ 


2 77 ^ 


(1-y )B 


d. 


( 1 + 1 .) 

4(1- 1)2) B gS 


4 = _ kJT^H (1 . 

2 ( 1 - V 2 ) B 


e_ = 


JT 


%2 


^ 32(1- i>2) "" 


(9 ^ 

B^ 


. = (g + 2_ ^ 9.5 

32(1- V^) B^ B^ 


(2.51) 
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^ I2(l+v) g2 


(1-V^) ISB^ SB^ 


. = rn + 1-^ 

'A " Q '' 1 O ' 

2(1- V^) B^ 


JT'^H 




2(1- v dB B" 


®6 


iJifL (1 + 1-) 

2(1- v^) B^ 


'7 


+ 1) 


2(l-y^)B B^ 


12(1+ V)H 


Q = 5 rr 

12(1+V)B,H 

neglecting in-plane inertia terms in eqs.(2,50)5 
combining these in a way similar to that in the case of 
the rectangular plate, and neglecting terms containing 
fourth and higher order deriva,tives with respect to 
time, the following equation is obtained : 

g^P"p2+g2P"P+g3P'‘+g4P(P»)^+g5(P‘)^+ggP^+g^P^+ggP = 0 

(2.52)* 

Comparing eq.(2.l6) with eq.(2.52), it is seen that qua- 
dratic nonlinearity is not present in the former. 

Its implication will be discussed in Chapter III. 
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'here 


§1 
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gg 

- 2(a^ 
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§3 

aib^ 

^2*^2 ^ 


- ± X ± s 7'’4 

age 8 - 12368 

6(ai '*■ + ^4^1 "533> 

2(a^ + h3_) (eg + e^3g + 653^ + + e^jg) 

(a^h3_ - aghg) (e^ ^ e^3, - 3^33) 

"■ ^2^2^ ^®2 + ^4^2 ^5^4 ^6^1 + ^70*3) 

(aib^ - aghg) (03 + e^jg + 6734) + eg(a2b3-a3b3L) 
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033.3 - C3di 
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CVI 

h 

Ogd^ - 

c^di 

“hi - 

‘^2*^2 

Q3'^2 - 

Cidg 

cidi - 

C2<i2 

^2 - 

^1*^4 


^ 1^1 “ 02^2 
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For small amplitude oscillations, t ems containing 
products of P and its derivatives are neglected in eq.(2.52}5 
giving , 

g3 P” + ggP = 0 (2.54) 


Corresponding non-dimensional period of flexural 
vibration is 


2 TT 


T 


^ yii/' rg 


(2.55) 


It can be seen that equations for the case of 
a rectangular plate derived in section 2.1 are obtained 
by putting c equal to 2 ero in the corresponding equations 
of the present section. 


2.4 THE case OF A SHALLOVf AFCH 

A shallow arch (see Fig .5) whose middle surface 
is given by ¥ = c Sin is considered here, c gives 
the initial maxd-mum rise in unstressed condition. For 
the arch to be shallox-/ c/a should not be greater than 
1/5. All the equations for this case can be derived 
as a special case by omitting variables connected with 
y direction and putting v equal to zero in the corres- 
ponding equations for the doubly curved shell given in 
the previous section*. 

* E 

This has to be done before ' g - ; - "-" — r 

for 0. 2(1+ V) 


is substituted 
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The stress resultants are 


X 


Eh (u^ ^ + 1/2 (w + 


C TT 


O jX 


X 


w ^ Cos 

• A. 


JT" x> 


_ Sh 
*'"x ~ 12 


X 


(2.56) 


5Gh 
= 6 ' 


( 0«..+ "W 2;) 


The boundary conditions along x = 0 and x = a 


are ; 


either 
“( = 0 


Uq = 0 

w = 0 


or 


= 0 
N = 0 

X 


(2.57) 


c jr 


<3^ 4-«,x 


'■x 


=0 


a 


For an arch hinged at Idth the ends and vibrating 
in the fundamental mode, rotation and displacements are 
assumed as, 

>s 


= A, (t) Cor 


*1 


u^= M(t)h Sin 
o a 

V = P(t) h Sin -JI-S 

a 


(2.58) 


Ultimately the second order equation, corresponding 
to eq.(2,52) of the previous section, comes out to be : 
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g^p»p2 + ggP^P + ggP” + g^PCP')^ + ggCP')^ + ggpS 


g?^ + §8^ " 0 


where 


=1 


3 4 p 

„jr h2 


3 4 0 

g2 = 2 * ^ ^ 


g = 2 5 STT^ 

3 {1+V)h2 12(1+V) 


g4 = 2 gj 

:z: rr 

5 

6 


+ 1 k2 JT^.H^ 

2 


= §2 


jr^H^ 


5 Ji 


4 


4(1+ ^t) 


o p 15 
5 H^k + — 


.6 




2 ^4„2 


) 


(2.59) 


(2.60) 


2 ' (1+V)H^ 

For small sjuplitude oscillations, terms containing 
prodncts of P and its derivatives are neglected, giving, 
§3^” + ggP = 0 (2.61) 

Corresponding non-dimensional period of linear 


vibration is 
2jr 


= 


■/Wi 


(2.62) 


Comparing eq.(2.4l) with eq.(2.59) it is seen that 
quadratic nonlinearity is not present in the former. 
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2.5 TtiE case of A SPHERICAL GAP 

Axi symmetric (torsionless) motion of a spherical 
cap has heen analysed in this section. The eq-uations of 
motion given by Ogibalor^" for the case of a shallow 
spherical shell (Fig. 6 ) are taken in this work after 
modifying them to suit the sign conventions used here. 
These force and moment equations were obtained by suitably 
integrating through the thickness, the equations of motion 
derived from the theory of elasticity. A spherical shell 
is considered to be shallow if ^ 4 (Ref .47). For such 
a shell, the actual values of forces and moments can be 
taken to be the projections of these in plan. 

The equations of motion are ; 

^rr “ ^Q<5 0Q p, 

^rr r r 5 (2.63.1) 

rr,r r ,r r 


Q 




, r p 


rr 


1 


- 5 (Tj.y + 0Q + Tae ) = Ph w 


rr 


(2.63.2) 


M Hh 

rr^r 


M - M 
rr 


se vPh2 i* 

■ - Q = “it" 0 


where u^ (r,z,t) = u^(r,t) + z, 0 (r 5 t) 


(2.63.3) 

(2.64) 


Q = k„Gh (w + 0) 
^ ,r 


(2.65) 
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In eq.(2.65) kg is a correction factor introduced to 
take care of distribution of shear stress through the 
thickness. For the case of a plate, its value has been 
derived by Mindlin to be by adjusting the lowest 

simple thickness -shear frequency of a homogeneous plate 
to its exact value. This value will be taken in all the 
comput at ions he re , 

The force-displacement and moment- curvature relations for 
the case of a spherical shell are taken as (Ref .62) : 


Bh 


•rr 


(l-v2) ” 


o o 


Eh ,0 ^ ^ 

Myy = I>(k^^ +y 

= D(k.,6,+i^ k^P 


( 2 . 66 ) 


’he strains and curvatures are taken as (Ref .59) 


j. iL 

rr r,r r 




- r + i 


krr ^ 


(2.67) 
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These are valid for moderately thick shells considered in 
this work; the transverse deflection is large (upto 2 or 
3 times the thickness of the shell), elongations and shears 
are small hut rotations are moderately large. Assuming 
the shell to he hinged along the periphery, w, 0, and u.^ 
are taken in the following form : 


w = A^Ct) h So (A r) 

0 = B^(t) (Xr) (2,68) 

u = C(t) (™ - r) 
a 


where , 

Sq ( Ar) = Jq( a r) - f I^( Xr) 
( A r) = J^( A r) + f I^( A r) 


(2.6S) 


These functions are assumed on the basis of linear 
classical solution of a circular plate. 

The boundary conditions, at r = a, for a hinged 
shell are ; 


w =0 

(2.70.1) 

= 0 

(2.70.2) 

li 

o 

(2.70.3) 


Boundary condition (2.70.3) is satisfied by the 
function selected for Up. The two constants X.a and f 
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are obtained by substituting the expressions of w and 0 
in eqs. (2.70,1) and (2.70.2). 

i'Jriting the equations of motion in terms of 
rotation and displacements, neglecting in-plane inertia, 
and applying Ga.terkin technique, the following equations 
are obtained ; 


^2^1 ^4;'^1^1 ®'5'^1 ■“ 0 ~ 0 (2,7l,l) 

’’hi + Vi + ”3^1° ■" ’’hi ” 5 ^? ”6®! ” 7 ° 

+ = 0 (2.71.2) 

b” + - CgA^ = 0 (2.71.3) 

xfnere , 


=<2(l+y).H + (l-v) H.lc„ j. X,. R* 


^ 1 

ag = - UnLiia . X, .R» 

^ A a 1 


ag = 2fi2 (X 2 + . X 3 ) 

(1- V )ks .H 

"■4 = ■■■ ■■ (X3 + 2.X2) 


as 




( A a)^ 


(X 3 + 2 .X 2 ) 



(1-v )K 
*- 

2 

(1- V )kc 


(X4 + X5) - 2 ( 1 + y ) .x^g, ( r *)2 


2.H. X 


(X. + x^) 


X + (1+ V) x,p + V Xq + 2 .x. 


^12 ■ " ^9 ■ “-^10 


X, 


11 


(l+y ).H.R* (X7 + '— + X,.) 


"IS- 


- (S.Xg + X^4) 


(1-- i> )kg.X^^.R* 

2H. ( X a)^ 


lliiLl .R*.(x^g + X^^) 


H 


( 1-^ ) kg.X^.^ 

2. X a 


. R 


- d-v'^) . X 


16 


" " (X.Q + X^^ 

(1-V ).Xt 


6k. 


^18 


Xg ) +- ^ 

^ (1+V)h2 


6. A a. k. 

(l+V )H 


ilJ KANPUR 

CENT i/j. ua«ARy 


'‘•v Ji 
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X. 


1 


So,. 


X. 


o 

a 


o 

a 


''^o 5 1* * j i*r *”ai*) d.3? 


X. 


^■^o,r) ^^“9-) dr 


X/ 


cl I S^^.S^ dr 


■'o*‘-o,rr 


o 

a 

„ f ^o*^Otr , 

a j — -2— {ir 


’^o* S^jrp-(2r«a)dr 


Xr- 


aj . S 
o 


o , rr 


dr 


X, 


8 


'^o * ^ r^ * S . dr 
° ^ 5 o 5 rr 


o 


X. 


'9 


(r-a) dr 


0 

a 

F 


X. 


10 


Sq. dr 


(2.72) 
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a/S,(s ) dr 


0 o,r 


®o-®o,r <2-? 5 to 


• Sq , r , 

a I — ^ dr 


3 r 


S-^ . (2r-a)dr 


1 f 2 


a Sn dr 


0 

a 

1 f 


■g S^.(r-a)dr 


-a (S^ _) dr 


=■ ®o,rrr- ^ o,t 


q g 

•^o-rr* o,r 


1 — 


o 



S4 


X21 


j2 


dr 


It does not seem to be possible to get one 
equation in though, however it is possible to obtain 
one equation in Combining the three equations (2.7l), 
neglecting derivatives of of the order higher than the 
second, the following equation is obtained s 

diBi + + dgB^ + + dgB^B![ + d B^" = 0 


Making use of eqs.(2.7l) and (2.73) 


(2.73) 


c.B. 


C2(d4+d5B;j_+d6B2) 


where 


1 


"^1*^1 a-,br7C-, 

“T™” + bo + 

°2 2 cg 


3-2^7 


a = Shh + a4>>70i , . ^ 


Cg “o~/ 2 ■ 

^2 

2 


+ 


h2l + ^ agbgcq 


(2.74) 


? 2 

_ a4b3cq a5b3Cq b^c 


dg ^ ^ 


1 




(2.76) 
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= 


. "*^9^1 

4 * ^ 


a b„ 
1 7 


= ,^^3^7 a^t>7 Sb^c^ bg aa^bgc^ a^bg 

Q -f ---«— *t- Q 4- -r 4- - 

Cg ^2 Cg °2 C 

= 2a-3^3^1 . %'^3 . 

^ -f 4- TT o 4" ■ - 


=1 


0'’ce is solved from eq.CB.VS), can be 
obtained from the cubic equation (2.74). 

For small amplitude oscillations, terms contain- 
ing products of Bj^and its derivatives in eq.(2.73) are 
neglected, giving, 

'^1®! ^4^1 = 0 (2.76) 


Corresponding period of linear oscillations is ; 

2 TT 

7 = - (2i77) 

'/iTT 

2.6 THE CASE OF A CIRCULAR PLATE 

The problem of a circular plate can be reduced 
from that of a spherical cap discussed in the last 
section. On putting R* to be zero, the final governing 
equation (2.73) reduces to : 

d^B^ + dgB^ + dgB” + d^B^B^ = 0 (2.78)'*’ 


4- 

Comparing eq.(2.78) with eq.(2.73), it is seen 
that quadratic nonlinearity is not present in 
eq.(2.78). 
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Where 


, _ ^ 
d_, - — — + 13 

-L Co ^ 


3 2 3 

^3^3^! . a4b3C3_ agbgC^ h^c-^ 

°i 4 =2 


h , '^9°! 


<=2 Cr 


(2.79) 


2a4b3C^ 3a3b3C4 a5b3 3b5C^ 

.2 3 — 3~ 




■r c^v 


are defined as in the 


the previous section. 


The period of linear oscillation is : 


2 rr 


(2.80) 



CHAPTER III 


RESULTS BID DISCUSSIONS 

In the previous chapter, a set of vibration 
problems has been reduced to that of solving second order 
nonlinear ordinary differential equations. In these 
equations non-dimensional time occurs as the independent 
variable. These equations have been solved by numerical 
integrations using the Range -Kutt a- Gill method. Compu- 
tations have been carried out on ISM 7044. The Poisson^ s 
ratio has been taken to be 0.3 in all the calculations. 

Table 1 contains the following results for 
square plates vibrating in the fundamental mode : 

i) periods given by the classical theory (Tq) 

ii) Mindlin's values* (Ref .63) for linear periods 
with the effects of transverse shear and 
rotatory inertia 

iii) linear periods (T^) with both the above 
effects included 

iv) nonlinear periods (Tpf) • 


Time periods have been computed from the equations 
given by Mindlin and shear coefficient has been 
taken to be sn-2/i2 as derived in Ref. 6l. 



this approach are a little higher. 
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Table 4 gives the values of periods of a beam 
of rectangular section for different thicknesses. The 
following results are given % 

i) periods given by the classical theory 

ii) periods for the linear case obtained from 

Timoshenko's equation^ taking shear coeffi- 

2 

cient equal to JT /i2 for the fundamental mode 
as derived in Ref. 9. 

iii) periods for the linear case (Tj^) as obtained 
in the present work 

iv) periods for the nonlinear case (T^J . 

i’4 

There is a close agreement between the linear periods 

obtained here and those obtained from Timoshenlco* s equations. 

Fig. 8 shoxxTs a plot of T /T as obtained in the present 

work. As in the case of plates the nonline arit 3 '' in this 

case also is of hardening t 7 /pe. For comparison the curve 

17 

given by Chu and Herrmann"" for thin beams, is also 
plotted*. The agreement between their results and the 
present results is good for thin beams (H^ = 30 or more). 
Tables 5 to 9 list the values of periods for the 
Doubly curved shell discussed in section 2.3 of Chapter 
II. The values are for different initial rise of the 
shell on square plan and for various thicknesses, 

* Perturbation solution given by Chu and Herrmann has 
been taken. 
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Values for various amplitude ratios are listed. Corres- 
ponding linear periods are also given in all the tables. 
Comparing Table 5 with Tables 6 through 9 it is seen that 
nonlinearity is of hardening type (similar to that of a 
flat plate) for a shell of slight initial rise (k = 0,2) 
whereas tbe nonlinearity is of softening tyye for shells 
of sufficient initial rise ,(k = 1, 2, 3 etc.). The 
elements for which nonlinearity is of hardening t 3 rpe are 
termed slightly curved elements and those for which 
nonlinearity is of softening type are termed shells. In 
the present case the transition from a slightly curved 
element to a shell takes place for k somewhere betv;een 
0,2 and 1,0, The phenomenon of change of nonlinearity 
from hardening to softening can be explained from the 
governing equation for plate vibration, eq, (2.16) and 
that for the shell vibration, eq, (2.52). The term 
which is primarily responsible for changing nonlinearity 
from hardening to softening is present fu the latter 
equation. This term is essentially associated with curva- 
ture terms and vanishes for flat elements. This has not 

"1 q 

been pointed out by Chauhan and Ashwell ^ , who studied 
vibration of thin Doubly Curved shells and flat plates. 



Figures 9 through 12 show the plot of the ampli- 
tude ratio versus period for shells of different initial 
rise . It is noticed that unlike in the case of a flat 
element the positive (outward) and negative (inward) 
amplitudes of vibration for curved elements are not the 
same. The negative amplitude is always larger. Dynamic 
snap-through is clearly seen from the figures. Before 
snap-through occurs the shells have nonlinearity of the 
softening type whereas after the snap-through the non- 
linearity changes to hardening type. Results for thin 
shallow shells of similar geometry are given in Ref. 3.9 
j,n the form of graphs. The trend of results is similar 
to that obtained in the present work. An accurate compa- 
rison with the present results is not possible because 
no tables of values are given. However, some limiting 
values have been compared and the comparison is good. 

For example, for a thin shell (say = 30) on square plan 
and with k = 1, the linear time period, T.;^- in the present 
work is 15.805. The corresponding p (the frequency parameter 
plotted on x-axis in Ref. 19) is 39.4ii*, which is nearly 

the value read on Fig. 6 of the said reference. For 
thicker shells, say Hj, = 20 and 10 the corresponding period 

* Period T in the present work is related to p as 

-r = 4 JT y3(l-V ^) 
p . H 
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as converted from this p will be 10.536 j 5.268 
which deviate from the present values of 10.560 ? 5.352 
respectively. The deviation is more for the thicker 
element. Similar comparison can be made for a shell on 
Square plan with k = 2 and k = 3 for linear periods and 
also for periods, at some other amplitudes. 

Tables 10 and 11 list the values of periods for 
different initial rise and various thicknesses of the 
Doubly Curved shell on a rectangular plan (B = 0.5). 

The trend is similar to that for a shell on square plan. 

Tables 12 through 15 list the values of periods 
for different initial rise and various thicknesses of 
shallow arch. Figures l3 through l6 show their plots. 

As expected the behaviour of an arch is sinular to that 
of the Doubly Curved shell discussed earlier. It is 
seen that the transition from a slightly curved beam 
to an arch takes place below initial rise equal to the 
thickness of the arch (k = l) . Dynamic snap-through is 
clearly seen from the figures. 

Values of the integrals through 
in section 2.5 dealing with spherical cap problem, are 
given in Table 16. Results for spherical caps are given 
in Tables 17 through 19 and those for a circular plate 
in Table 20. Figures 17 and 18 show plots of period 


Xg^, defined 
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versus amplitude ratio for spherical caps of R* = 0.1 
and 0.2 respectively. Shallow spherical caps behave in 
a way similar to other shallow elements discussed earlier. 
The fact that the positive negative ampl J.tudes of vibra- 
tion are different j.n spherical caps as in other shallow 
elements has; not been poj.nted out by Grossman et. al.^"^ 
who have given results for thin spherical, caps. Also, 
thev have not discussed the phenomenon of dynamic anap- 
through. Figure 19 shows a plot of Ti/T^p versus R* 
for spherical- caps (Tpp stands for linear period of a 
circular plate). In this figure a curve due to Grossman 
Gt.al. is also given ^ For =30, the curve based on 
present results almost completely coincides with their 
curve. Figure 20 shows a plot of Tr/T^ versus Ap for a 
circular plate. Nonlinearity is of hardening t;';pe and 
behavi-our of a circular plate is similar to that of a 
rectangular plate. The transition from a slightly 
curved plate to a shallow shell occurs at some value of 
R* below 0.1. 


+ 


The graph has been plotted here after suitably 
'transforming the x and y variables of the curve given 
by Grossman et.al. 



CniPTEE lY 


GO'NCLUSIOKS 

i) It is known, that in the case of small amplitude 
flexural vj.b rations , transverse shear and rotatory inertia 
decrease the frequency. This trend continues to be so 
even for large amplitude vibrations. This is seen for 
all the structural elements analysed in this work namely, 
straight and curved beams, rectangular and circular 
plates, a Doubly Curved shallow shell and a spherical cap. 

ii) For any flat element like beams and plates, the 
nonlinearity is of hardening type (that is, the period 
decreases as ariplitude increases) whereas for shallow 
elements like shells and arches, the nonlinearity is 

of softening type (the period increases with amplitude). 

A distinction can be made between a slightly curved 
element and a shallow element by determining the value 
of the initial rise beyond which the type of nonlinearity 
changes from hardening to softening. This happens because 
of the presence of quadratic nonlinearity tenn in the 
final equation for curved elements. 



65 


iii) In a shallow element, when some initial amplitude 
is given, it vibrates about the initial equilibrium 
position and the positive amplitude (in the direction of 
the initial rise) is always smaller than the negative 
ampliti-ide. The nonlinearity is of softening type, 

VJhen a critical value of initial amplitude is 
given to the shallow element, it snaps through to the 
other equilibrium position and stays snapped inside out. 

At thi.s stage, the frequency of transverse vibration 
becomes zero (that, is, the time period becomes infinite). 
When an amplitude of magnitude greater than the 
oriti,cal value is given, it passes through both the 
equilibrium positions during its oscillations and the 
negative amplitude is very large. The nonlinearity 
changes from softening to hardening type. It is assumed 
that the material remains elastic throughout. 

iv) The period of vibration decreases as Initial rise 
of a ciirvBcl olenorrt incroasos* 

v) Of particular interest is the feature that the 
nonlinear effect counteracts the curvature effect. It 
has already been observed that the period of vibration 

of a curved element is less as compared to that of a flat . 


element. 
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The percentage decrease is less if the nonlinear theory 
rather than the linear theory is employed.* 

lor example, a straight beam and a shallow arch 
of initial rise equal to four times the thickness (that 
is, k = 4) are considered. Hatio of length of span to 
thickness in both the cases is taken to be 20 (that is, 

Hy = 20) and the nonlinear periods of vibration for an 
amplitude ratio of 1.0 are considered c The period based 
on nonlinear theory decreases from 24.77 to 4,80 whereas 
that based on the linear theory decreases from 44.29 to 
4.49. The percentage decrease in the period given by 
the nonlinear theory is much less than that in the period 

m 

given by the linear theory. 


The comparison of periods is made before the dynamic 
snap- through occurs. 



CHAPTER V 


RECOMENDATIONS FOR FURTHER WORK 

i) In the present work the transverse deformation has 
been taken as ; 

w(x,y,z,t) = ¥(x,y,t). 

Thus the effect of normal strain is neglected which is 
justified because the fundamental frequency of flexural 
vibration is much less than the frequency of thickness- 
stretch mode of vibration. However, the effect of 
normal strain can be included by taking w in the form : 
w(x5y,z,t) = WQ(x,y,t) + z.f (x,y,t) 

ii) Sine functions have been assumed in this work for 
the initial unstressed shape of the Doubly Curved shallow 
shell and also for the shallow arch. Different geome- 
tries can however, be analysed. The effects of imper- 
fections in the initial geometry on phenomena like 
dynamic buckling etc. can also be studied. 

iii) In this work, elements with hinged ends have been 
analysed. Extension can be made to other boundary 
cond.ltions. 
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iv) The coupling of modes between the fundamental and 
higher ones has been neglected in the present work. The 
mode shape can be refined by taking additional teims. 

For example, for a hinged beam vibrating in the fundamental 
mode, the tronsverse deflection can be considered as ; 

r n-y: 3 xrx 1 

w = P(t) Sin + P-, Sin + 1 

L- < a ^ 

and similarly for other components of displacement and 

rotation. The additional unknowns p^, etc. can be 

solved from the additional equations obtained from 
variation; that is, by equating the coefficients of 
^P]_? zero. 

v) Free vibration problems have been analysed here, 
however, forced vibration problems can also be studied 
in a similar way. 

Vi) Extension of the present work to deep shells is 
feasible. In that case, the displaceraent and rotation 
functions should be selected judiciously and any inte- 
gration to be carried over the surface of the shell 
should not be carried over the projected plan. 
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vii) If the amplitude of vibration is not small^ some 
materials may exhibit nonlinear elastic properties. In 
such a case, the present work should be extended to 
include both geometric and material nonlinearities. 

viii) In this work, the structural material has been 
assumed to be isotropic. However, anisotropic materials 
are used widely in the modem age. Extension of this 
work for such materials is possible by taking suitable 
constitutive equations. 
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table 5 

PERIODS FOR A D.C, (DOUBLY CURVED) SHSLL(B=l|kr: ,0 .2) 
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In this and subsequent tables, figures in brackets given below the values of 
T|^ indicate corresponding negative (inv/ard) amplitude ratio. 
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TABLE 10 

PERIODS FOR A D.C. SHELL (B^O.S-, k=l) 
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TABLE 11 

PERIODS FOR A D.C. SHELL (B=0,5;k=4) 
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TABLE 12 

PERIODS FOR A SHALLOV/ ARCH (k=l) 
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^PERIODS FOR A SHALLOW ARCH (k=2) 
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PERIODS FOR A SHALLOW ARCH (k=3) 
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PERIODS FOR A SHALLOW ARCH (k=4) 
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TABLE 16 


SOME BITEGRALS IM SPHERICAL CAP PROBLEM 


Integral 

Value 

Integral 

Value 

h 

0.16959 

X 12 

0.46258 

^2 

-0.20747 

^13 

-1.07608 

^3 

-0.35158 

Xi4 

-0.97503 

X 4 

-1.11781 

45 

-0.16959 

^5 

-1.38822 

^16 

0.45219 


0.51098 

Xi7 

-0.38593 


-0.93106 

^18 

1.11810 

^8 

-0.59501 

Xi9 

-3.36817 

X 9 

0.81439 

^20 

3.29342 

^10 

-0.46282 

X 2 I 

4.65301 

^11 

0.46553 




( A a and f have been calculated to he 2.2174 and 0,0379 
respectively. These values have been used in evaluating 
the above integrals). 


table 17. 

PERIODS FOR A SPHERICAL GAP (R*=0.l) 
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TABLE 18 

PERIODS FOR A SPHERICAL CAP (R*=0,2) 
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TABLE IS 

PERIODS FOR A SPHERICAL CAP (R* =0.3) 
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PERIODS FOR A CIRCULAR PLAT 
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FIG. l_GEOMETRY OF A RECTANGULAR PLATE 
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FIG. 5_ GEOMETRY OF AN ARCH ( UNDEFORMED ) 



FiG.6-(a) GEOMETRY OF SPHERICAL CAP ( UNDEFORMED ) 
(b) STRESS RESULTANTS IN SHELL ELEMENT 
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FIG.8_ T,^/Tc versus Ar FOR A BEAM 
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FIG. 17- PERIOD VERSUS AMPLITUDE RATIO FOR A* SPHERICAL 
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